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I. Introductory, 

In his important observations on the lunar deflection of gravity, Dr. 0. Hecker 
has pointed out that the force acting on the pendulum at Potsdam is a larger 
fraction of the moon's force when it acts towards the east or west than when it 
acts towards the north or south. Various explanations of this anomaly have been 
proposed ; among them one, suggested by Prof. A. E. BL Love, # is that a possible 
cause may perhaps be found in the attraction of the tide wave in the North Atlantic. 
Recently Prof. A. A. MichelsonI* has found a similar result in his arduous task, 
successfully achieved, of obtaining a continuous record of the lunar perturbation of 
a very long water-level at Chicago. Prof. Sir J. Larmor kindly suggested to me 
a query, whether the excess-pressure of the tide in the North Atlantic would affect 
much the measurement of water-level at Chicago, owing to the elastic depression of 
the earth's surface that it would produce. 

In the present paper the calculation is undertaken in order to ascertain to what 
extent the consideration of tilting of the ground is important for the explanation 
of this geodynamical discrepancy, on the assumption that the earth consists of elastic 
matter of uniform density and of uniform invariable elastic constants, and that the 
North Atlantic may be represented roughly as a circular basin and that the tide in it 
is uniform or else elliptic. The curvature of the earth is neglected, which is, of course, 
admissible for a first estimation. 

II. General Formulae. 

The solution of the equation of equilibrium concerning an elastic body with plane 
boundary, otherwise extending to infinity, subject to certain boundary conditions, 

* 'Some Problems of Geodyimmics/ Cambridge (1911), p. 88. 
t ' The Astrophysical Journal/ vol. 39 (1914), p. 105. 
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has been obtained by Lame and Clapeyron* and by J. BoussiNESQt ; the former 
by making use of Fourier's theorem, and the latter by introducing several kinds 
of potential function. The solution for the case where the boundary condition is a 
normal pressure distributed symmetrically round a point on the surface and the body 
is free from bodily force is given by Prof, H. Lamb J in a form of definite integral. 
I have tried to solve the problem, independently of Lamb, and obtained the solution 
answering to any distribution of a normal pressure, from which the solution of Lamb 
can be derived as a special case, § 

For the present purpose let us confine our attention to the simplest case in which 
the boundary is subjected to a symmetrically distributed normal pressure. Take the 
centre of the area, in undisturbed state, on which the given normal pressure is 
distributed, as the origin of the co-ordinates, the inward normal. to the undisturbed 
surface as the z-axis, and denote the distance of any point from the 2-axis by r. 
Then corresponding to the normal pressure 
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fix), . . . . (1) 



given at the surface z = 0, we shall have the following expressions for the displace- 
ment components :— 
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in which Z(k) stands for 



Z(k) = k\ f(x)J Q (kx)xdx, ........ (4) 

J 



X, jul being Lame's elastic constants, J (x), J x (x) Bessei/s functions of zeroth and 
first order, and a the radius of loaded area. 

Now, suppose the normal pressure (l) to arise entirely from the tidal loading in 
the North Atlantic basin, then the tidal loading would come into play to affect the 
water-level measurement at Chicago in two accounts : the one is the direct attraction 
of the material load, the other the deformation of the ground by the pressure 

* 'Crelle's Journal, 5 vol. 7, p. 400 (1881), 

f ' Application des Potentiels / Paris (1885), 

t i Lond. Math. Soc. Proc.,' vol. 34, p. 276 (1902). 

§ The solution, with numerous examples, will be published shortly elsewhere. 
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produced by that load. Dr. C. Chree,* and afterwards more completely Prof. H. 
NAGAOKA,t find: a formula, by using the formula obtained by Boussinesq, to calculate 
the deviation of the direction of gravity due to the attraction of a material load on 
the surface of the earth. 

The same result can be attained, of course, from our solution. The expression for 
the vertical displacement at a point on the surface can be transformed into 

(u z ) {) = — • —. -t — z-^jJ- x dx d</> 



2?r 2 (A-f/x) /a 
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o J 



K' 



by making use of Neumann's addition theorem for Bessel's functions, where R' 

stands for 

B/ = ^/(r 2 — 2rx cos <p + x 2 ). 

On the other hand, if we denote the attraction constant by y, and gravity, prior to 
the application of the load, by g, then the gravitation-potential at a point on the 
unloaded surface due to the loading can be expressed by 

* 

V == y ■ - • ■- L % T l X dx d(j>, 
J o * o g \\> 

provided the height of the loading material is negligibly small compared with the 
distance of the point under consideration from any point in the loaded area. 
Comparing the above two expressions, we have 

Vu " g~ A + 2 M [ z)() ' 

Thus the direction of gravity becomes in consequence of the attraction of the 
loading material inclined to the vertical at the angle \/r which will be determined by 

while its tilting effect is expressed by 

***-&). <•> 

The total effect of the loading will thus be 

^H^-^XIi (7) 



in a close approximation. 



* 'Phil. Mag./ V., 43, p. 177 (1897). 

t Tdkyd, <Sug. Buts. Kizi,' VI., p. 208 (1912). 
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III. Uniform Loading. 

Let a be the radius of the loaded circle, h the height of the material loading, which 
is supposed to be uniform, and p its density, then we have to put 



f(r) = —ghp for r < a, 



= 



?? 



r > a. , 



On this supposition, we get 



therefore 
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The other components of the displacement and those of stress can be expressed in 
similar forms. But it is unnecessary to write them down here as they are out of our 
present purpose. 

The integrals required here cannot be evaluated in a very simple way. Some of 
them are closely connected to the magnetic potential due to a circular current, or to 
the velocity-potential and stream function of a circular vortex, and have been 
discussed by various authors. In his paper on the inductance of circular coils, # 
Prof. Nagaoka has devised a comparatively simple method which may be applied to 
evaluate all the integrals needed for the calculation of the displacement and stress in 
the present problems. Let us follow his method and describe it here briefly. 

Put 

R — v / (a 2 -\-r 2 —2ar cos 0), 

then by Neumann's addition theorem for Bessel's function we have 
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Phil. Mag.,' VI. , 6 (1903), p. 19. 
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Making use of these formulae we obtain 
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To find these integrals, put 
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so that 



03 < «2 < 6 1? 

ei + e 2 + e 3 = 0, 
and change the integration variable from (9 to 5 by 
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Put again 



s = JP (w), 



then, since 8 or P (u) is real and lies between e 3 and e 2 , s = e 3 and 5 == e 2 correspond to 
w = w 3 and u = w 2 respectively, if we take 9 f (u) to be positive # ; where w x and a> 3 
denote the real and imaginary half-period and oo 2 = ft>i + a> 3 . Thus 



,7r cos e de 
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For the evaluation of the other integrals, write 



9 (v) = 



2ara 



. . . * I Xo) 



• • a 



• • 



(14) 



* If we assume P (£4) to be negative, then s == e% and s =e% correspond to u ~ &% and u » 2<*)i + w 3 
respectively. But the same result will, as a matter of course, be obtained after integration. 



40 



MR. K. TERAZAWA ON PERIODIC DISTURBANCE OF LEVEL 



then we shall have 
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Now we have 
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The term miri enters because of the many-valued property of a logarithm. 
actual value of m and 9 f (v) will be determined by the following consideration. 
From the definition of P (v) and e l9 e 2 , e 3 , it follows Immediately that 
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accordingly 



e 2 <f>(v) <e } 



The last inequality shows that the value of v must be one of the following : — 



(i.) v = (2n+l) w 1 -{■ (2n f + d) w s 
(ii.) v = (2^ + l)a> 1 + (2^ / + 2-0)a)j 






>• 



(18) 



where n and n' denote any integers, positive or negative, or zero, and 6 a positive 
number less than unity. To determine the value of m in the formulae (16) (17) for 
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the value of v given in (i.) of (18), observe that the integral on the left-hand side of 
(15) and the function v^ — co^fo) change their values continuously as 6 varies from 
to 1, while m remains unchanged in this variation. In the limit as 6 -> 0, the value 

of the integral is nil and 

2^ — 2w 1 £(v) = 2n f 7ri, 

and therefore we have 

(i.) m ==— -n'. 

Similarly for the value of v given in (ii.) of (18), proceeding to the limit 6 ■> 0, we find 

(ii.) m = — (n f +l). 
The value of p' (v) will be obtained from 



P' a (v) = 4 [P (v)-ej . [P (v)-ej . [P fa)-e 8 ] 



In the present case we have 
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It may easily be shown that the value of ¥ (v) is a positive imaginary quantity for 
the value of v given in (i.) of (18) and a negative imaginary quantity for (ii.). 
Therefore we have to take 
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Lastly, one more integral which can be evaluated without the knowledge of elliptic 
functions is 
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Substituting these values of integrals in (9) ? we have 
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in which the simplest value of v, i.e., n = 0, n' = iij, (i.) of (18) is taken as the 
representative. 

At the surface it reduces simply to 
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Similarly for the value of ( ^~ ) we have 
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For practical purposes of calculation, it will be very convenient to transform the 
expressions into Jacobi' s q series ; q being defined by 
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then g can be computated from 
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At a distance from the edge of the loaded circle, I is a fairly small quantity and 
consequently the terms after second or third may be dispensed with. The q series of 
the functions needed here are as follows : — ■ 

7r 2 fl °° q 2n 

Wl = f(f)V(0), 

x _ i ry" ] (o)_y,(on. 

■ ***-*- 8^1^(0) &,(o)J' 
S- 2 (0) = 2q i (l + q 2 + q (i + q w +...), 

$"\(o) = -2ttV (i-3Y+5Y-7V 2 +-)> 

^(0) = 2xg l (l-3g 2 +5g 6 -7g 12 +...), 

At the point near the edge of the loaded circle, the above expansions cease to be 
applicable. For this case, our object will be accomplished by using the quantity q x , 
instead of q, which is defined by 
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Making use of the transformation formula of Theta-functions it will be easily 

shown that 

^0(0 |t) 9 . , 8^(01x0 
9o (0|t) ^(OItx) 

^(01t) „q, , 3 &" 8 (0 ] T t ) 
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an 
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consequently we have 
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The g x series for the functions required here are 
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In the accompanying diagram the courses of (u z ) Q and its slope are exhibited as 
functions of the distance of the point of observation from the centre of the loaded 




f— . 5> t 



Eig* 1. 

circle, with the proper scale. The lower curve thus represents the vertical displace- 
ment and the upper one the corresponding slope outside of the loaded area. 

If we liken the North Atlantic to a circular basin of a large radius and determine 
the relative position of Potsdam or Chicago referred to the centre of it, the attraction 
effect and the tilting effect of the periodic filling and emptying of tide, which might 
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assist in producing the extra east-west force, in observations of lunar disturbance of 
gravity, may be computed. If we suppose the place of the observation to be not very 
near to the circular basin, the effect, as we see from the above diagram, is of course 
small, but it increases rapidly as the edge is approached. 

For example, suppose the radius of the North Atlantic basin to be 2000 km., 
the position of Chicago to be 3000 km. from the centre, and the level of the water 
in this area to be raised one meter, then 

— = 1-5, g^ = 0'00255, 

a 

and 

aa 2 (k»r^) = 0*8639. 

Further assume that the density of sea water is 1, and in c.g.s. units, 

y = 6*65 xl0~ 8 , g = 980, 
X + 2/x 3 
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M = 6xl0 n , 



then we shall have 



2(X + m) 

V'= l'17xl0- 8 = 0"'0024, 
= 3'37 + 10- 8 = 0"'0069. 

Accordingly the total effect amounts to 



x/ r + = 4'54xl0- y - 0'"009. 

It will be noticed that the effect of tilting is about three times as great as that of 
the attraction, so far as the material constants are assumed as above. According to 
Lord Kelvin, # who initiated these investigations, the direct lunar effect on the 
deviation of a plumbline is a maximum when the moon is at the altitude 45 degrees and 
amounts to 7/, 017 nearly. The total effect of a tide of amplitude one metre (which 
is possibly two or three times the actual amount) found here is not small enough to 
be neglected compared with the direct effect of the moon, the former is nearly half 
the latter. As the tilting effect and the attraction effect of the tide wave are 
directly proportional to the height of the tide, the total effect oscillates in time in 
accordance with the law which the tide obeys. There is, in general, a difference in 
phase between the lunar effect and the tidal effect, which is worthy of closer 
investigation, t But we must bear in mind that the calculation adopted here is 
nothing but a rough estimation of order of magnitude, since the North Atlantic is far 
from circular, the tidal loading in it is never uniform. 

* 'Natural Philosophy,' Part II., p. 383. 

t Fortunately, the phase difference of both effects in Michelson's experiment may be neglected in a 
rough estimation, owing to the relative position of Chicago and the centre of the North Atlantic. 

H 2 
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Nevertheless the above analysis shows that the tidal effect on the water-level 
measurement, even at a point as far from the coast as Chicago, plays an important 
role and cannot be regarded as a small correction. 



IV. Elliptic Loading. 

Next, let us suppose that the tide in the North 
surface is given by the equation 

2 2 

2 a 2 . 



Atlantic is not uniform, but its 



viz., the excess pressure on the bottom due to the tide diminishes on approaching the 
coast so as to amount to 
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In this case the function Z(h) becomes 
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Therefore we have 
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The integrals contained in the above can be obtained by expanding the trigono- 
metrical functions into power series of h and making use of the formula 
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where 
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vV+z 2 ) 



These series converge for \/(z 2j rT 2 ) > a, and are applicable in this region. 
At the surface we have to put z = and v = 0. Since 
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for r ~ a. 

To find the expression for the vertical displacement within the loaded circle we 
proceed as follows : — ■ 

Making use of the power series of Bessei/s function, we have 
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The evaluation" of the function Q m (x) can be undertaken by the aid of the formula? 
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A little calculation will give us 
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and, in general, 



Q m (x) = (-1) 
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Thus the integral on the left-hand side of (33) can be expanded in an ascending 
power series of r/a which probably converges for limited values of r if the value of z 
is fixed. This series and that found in (30) have a common region in which they 
are both convergent and therefore they must be congruent to each other in that 
region. On the proof of this proposition we shall not enter, but we shall find the 
region of convergency of the latter series at the boundary. 



z 



Expand Q 2n ( — ) into a power series of z/a, supposing zja to be sufficiently small, 



z 



then the first term of it will be ( — l) n 2 (n — l) (2n — 2)\ — . If we retain only the 

a 

terms which contain the first power of zja in the series of (33), its general term will 
then be 

2(n-l) (2n-~2) \ (v\ 2n z_ 
(n\) 2 2 2n W a 

The series which has this expression as its general term converges obviously for 
the values of r smaller than a. Thus the expansion (33) applies for r <<x when z is 
an infinitesimal. 
Since, for z ~ 0, 



we have 



for r ~ a. 
Hence 



G (0) = J, Q a (0) = |, Q 4 (0) = 0, Q 6 (0) = ! 



>oo 



sin ha — ha cos ha 



Fa s 



3 



J (hr) dh 



7T i . i / r 



4a 



1- 



2 



\aj J 



/ \ = 7rabgp(\ + 2fi) l, x .( r \ 



(34) 



(35) 



for r = a. 

Quite similar arguments may be employed to find the expression for the tilting. 
We shall have 



• CO 



sin ha — ha cos ha 
¥a % 



J± (hr) dh = 



rv* 



r f • 1 a 
- snr - 

I r 



2d' 



a 

r 



1 



a 

r 



2\r> 



"\ 



j' 



r 



/~4 



r 



• (36) 



/i'» 



7T7 

4a' 



a. 



j 
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Consequently 



on, 

dr 



bgp(\ + 2jui) r f . ■ ,a a,.* 



4 / u(X + /x) a { v 



r 



9\ - 1 



■ • (37) 



for r 



€&• 



The general march of the function (u z ) and 



3-r/o 



are shown in the attached 



diagram (fig. 2) in a suitable scale. The lower curve represents the vertical 
displacement and the upper one the slope of the ground. The course of the curves is 
very similar to that of fig. 1, as we would expect, except at the points very near to 
the edge of the loaded area, where, as there is no discontinuity in the load in this 
case, the slope remains finite. 




> r 



Fig. 2. 



Let us take an example, with the same assumption regarding the various constants 
and the position of the point of observation as in the former example, and suppose 
that the total amount of the load is the same as before, i.e., the mean height of the 
tide is one metre. Then we shall have 



an 



d 



xj, = 112xl0~ 8 = ;/, 0023, 
£= 3'21xl0~ 8 = 0"-0066, 



^ + 9 S = 4;33 x 10- 8 = 0" # 009, 



nearly the same as the result in the former example. 
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If we suppose the place of the observation to be nearly at the edge of the loaded 

area, then 

f = 5'0xl0- 8 = ;/, 01, 

<j> = 14*4 + 10- 8 - 0"'03, 

and 

^ + ^> = o"'04, 

amounting to more than the direct lunar effect itself, notwithstanding that the load 
falls off towards the coast. 



Mjaaenua inmuwi 



